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1) , In this work we compute the leading logarithmic corrections to the 

6 — > S7 decay in a dimensional scheme which does not require any 
definition of the 75 matrix. The scheme does not exhibit unconsis- 
^ ■ tencies and it is therefore a viable alternative to the t'Hooft Veltman 

H , scheme, particularly in view of the next-to-leading computation. We 

confirm the recent results of Ciuchini et al. 
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I. INTRODUCTION 

The h —^ s'-f decay has received considerable interest in the last years, thanks 
to its sensitivity to physics beyond the Standard Model. 

On the experimental side, the limit for the inclusive decay based on the first 
observation of the B -^ K*'~f transition by the CLEO collaboration [1,2], are com- 
patible with the SM predictions, taking into account leading QCD corrections. 

On the theoretical side, after the first computations of the decay ampli- 
tudes [3-7] and the confirmation of the results of [3] given by us in [8], several 
improvements have occurred. 

In [9,10] the analysis has been extended to include all the relevant opera- 
tors that mix under QCD corrections and contribute to the leading logarithmic 
corrections. 

In [11,12] the long-standing problem of the differences between the results 
obtained in [3,8] and in [4,5] has been solved, by showing that the anomalous 
dimension matrix is scheme dependent even at the leading order. This effect is 
compensated by the matrix elements of 4 fermion operators, and the physical 
results are scheme independent. A similar analysis limited to a subset of the 
relevant operators is also presented in [13]. 

In [14] it has been performed an analysis of all the sources of uncertainty in the 
leading logarithmic order computation (LLO) and it has been shown that the in- 
clusion of the next to leading order corrections (NLO) can reduce the uncertainty 
from 25% to less than 10%. 

In view of these considerations we think that it is time to undertake a NLO 
computation and in this paper we reconsider the LLO computation with a tech- 
nique that we think particularly apt to this purpose. As a byproduct, we confirm 
the recent results of Ciuchini et al. [11,12]. 

Dimensional methods differ in the treatment of the 75 matrix: to our knowl- 
edge, the only method that does not exhibit any inconsistency is the t'Hooft Velt- 
man scheme (HV) [15]. Unfortunately, this method is very cumbersome and in 
literature simplified versions have been used, like naive dimensional regularization 
(NDR) or dimensional reduction (DRED) [16]. It has been shown ( [17,18,11,12] 
and references therein) that NDR, DRED, HV methods give the same physical 
results in a large class of radiative-correction computations (up to two loops) 
where no trace ambiguities arise; on the other hand it is well known that inher- 
ent pathologies of NDR and DRED schemes are present at three-loops level [19]. 
We are therefore interested in a practical approach to the evaluation of radiative 
corrections, as proposed in [20] , avoiding the extension of the 75 matrix to rf 7^ 4. 
This method appears to be free of ambiguities, and in our opinion presents less 
complications than the HV scheme, especially when performing computations by 
symbolic manipulation programs. 

The plan of the letter is as follows: in section (II) we describe the scheme we 
have used, recalling the results in [20], then in section (HI) we show how we have 



applied the method to the b -^ 57 process. 



II. STRATEGY 

We aim to evaluate the amplitude for a decay process induced by the exchange 
of virtual particles, such as the flavor changing neutral current processes. 

The integration of the heavy degrees of freedom (in our case, W boson and t 
quark) leaves an effective hamiltonian, whose matrix elements in the background 
of the "light" theory determine the decay amplitude at the lowest order in the 
heavy-mass expansion 

(/ 1 ^)fuii =^j:c, (m, ^) (/ io,i ^),^,, (^) + o (^) . (1) 

The coefficients of the operators in the effective hamiltonian are determined 
matching the two sides of (1) in perturbation theory at a scale of the order 
of the heavy masses. At /i ~ M all the QCD large logarithms originate in the 
operator matrix elements of the effective theory and the improved perturbation 
theory allows to resum them sistematically. The RG invariance of the l.h.s. of (1) 
is exploited, together with the calculable /x dependence of the matrix elements, to 
evolve down the r.h.s. at values of fi comparable with the scales of the external 
states; this amounts to "transfer" the large logarithms from the matrix elements 
to the coefficients. 

The key element is therefore the evaluation of the anomalous dimension matrix 
of the operators, which determines their /i evolution, together with the "match- 
ing" conditions at the M scale and the matrix elements at the lower end of 
the evolution. This procedure is in general dependent on the regularization and 
renormalization scheme chosen for the operators. 

Now what is important for our discussion is the freedom to define the d di- 
mensional extension of the operators of the effective hamiltonian: as in [20] we 
will use this freedom to simplify the computation. 

Let us first observe that the general structure of the effective hamiltonian 
for the 6 — >• S7 decay contains chiral projectors Pl/r = 1/2 (1 =F 75); therefore a 
definition of the 75 matrix is required in d dimensions. We have 

'Hl, = J:C^R^ + J:C^'iL®L)^ + J:cl^''{L®Rh , (2) 

i i i 

where Ri stand for the magnetic momentum operators, while {L L)- , (L (8> -R)j 
are the current-current operators: the field content is immaterial at this level. 
Consider on the other hand the effective hamiltonian obtained with an exchange 

75 -* -75 

nl, = E C^L[ + E Ct' (R ® i?): + E C^"" (R ® L)[ (3) 



^s/a = 7T (^cff =t ^cff) • (5) 



where we intend that the operator L'^ is equal to Ri apart the chiral projector, 
and analogously for the other operators. 

Since QCD does not know the sign of 75, the two hamiltonians will receive 
the same perturbative corrections, that is 

VN[L ® Li] + 7f/^Ar[i?,-] + 7f/^^Ar[(L ® L)^^ + 7f/^«Ar[(L ® i?)^.] = 
VN[{R ® R)[] + 7f/'^iV[L;.] + 7f/^^Ar[(i? ® i?);.] + 7f/^^Ar[(i? ® L)'^^ = 

(4) 

assuming the same renormalization scheme is used. The A^ [O] symbol stands for 
the operator O renormalized by minimal subtraction. 

Let us now take the linear combinations which are symmetric and anti- 
symmetric under the 75 sign flip 

1 
2 

The original hamiltonian, which determines the physical amplitude, is a combi- 
nation of the symmetric and anti-symmetric hamiltonian; the key point is that 
in order to obtain the scaling properties of Ti^g it is sufficient to know the QCD 
evolution of Tig. In other words it is possible to choose the same renormaliza- 
tion scheme for Tig and ?ia, and therefore determine, for instance, the evolution 
of an operator having chiral structure {V — A) ® (V ^ A) from the evolution of 
{y ®V)±{A® A). In the case at hand we can choose a scheme for the symmetric 
part which is as simple as possible, and forget the anti-symmetric part. 

Let us focus on the current-current operators; the extension of the sym- 
metrized hamiltonian Tig to d dimensions can be conveniently done reexpressing 
tensor products of 7 matrices as elements of the Clifford d dimensional algebra, 
for instance 

- [{V -A)®{V-A) + {V + A)®{V + A)] = 7^ ® 7^ + -7^,^ ® 7^,^ . (6) 
More generally a complete basis in d dimensions is provided by the operators 

(7(")®7(")) = Y. 7m,/.2,...M. ®7/.i,M2,...M» 



(- 



7Ml.M2,---/in ~ I 2^ \ ^) 7mi7m2 • • • 7/i„ • \') 



This extension to d-dimensions does not require any definition of a rf dimen- 
sional 75 and therefore it avoids the complications due to the splitting between 
e-dimensional space and 4-dimensional space. It also allows to treat relevant 
operators having n < 4 and evanescent operators having n > 5 at the same way. 

Many tecnical details will be explained in a more extended work [21] where 
a detailed account of the computation and the full set of formulas needed in 
addition to the ones in [20] to implement the method will be given as a reference 
for future works. 



III. THE LEADING LOG CORRECTIONS TO THE B ^ S-/ 

HAMILTONIAN 

Let us write down the effective on-shell hamiltonian for tlie b -^ 57, b -^ sg, 
already used in [9,11], in euclidean notation: 

Ql = {SaC^)y_^{c^ba)y_y^ 
Q2 = (SaCQ,)^_^(C/36/3)^^_^ 

Q3,5 = {Saba)v^A J2 {mp)v±A 
q=u,d,s,c,b 

Qi,& = (Sa&/3)y_^ Yl {Ql3qa)v±A 
q=u,d,s,c,b 

^ i-ieQd) _ „ , 

(47r) 

Qs = ———T''^bScr^j.u(V+A)G^,yb. (8) 

(4vr) 

Following the steps outlined in the preceding section, we consider the on-shell 
basis, in d dimensions 

0(1, „) = ^ (s«7^"^c^) (cpl^'^K) 

0(2, n) = ^ (s.7^"^C„) (c/37^"^&/3) 

0(3, n) = ^ (sal^^'K) E (^/57^"^g/5) 

q=u,d,s,c,b 

O5 = ,, ,2 rUbSa^^F^^b 
(4vr) 

Oe = , ^ ,2 rribSa^^G^^b . (9) 

(47r) 

The symmetrized operators corresponding to the ones listed in (8) are defined as 

Ql = 0(1,1) + 0(1,3) 

Q2 = 0(2,1) + 0(2,3) 

Q'3/5 = 0(3,1) ±0(3,3) 

Ql/e = 0(4,1) ±0(4,3) 

Qr = 05 

Qs = Oe. (10) 



The first step is to determine the anomalous dimension matrix of the operators 
hsted in (9), by using algebraic identities for arbitrary n [22,23,20]. 

The next step is to determine the Renormalization Group evolution in d di- 
mensional space: 



VN 



On 



(i, n) 



+ E7C 



i, n), (j,m) 



N 



J,m 



VN[0,,] +J2lnmN[0, 



0(, 




k 







;ii) 



where 



V 

Pqcd 



d ^ d d 

fJ'-K- + PQCD^:^ + Inn, rUb 

2e-g, 



dgs ' "'"""""dmb 
as 1 1 Ca — '^iT'f 



:2n) 



(47r) 
3Cp 



(12) 



here and in the following Ca, Cp are the Casimir invariants of SU{Nc), Ca = 
Nc, Cp = {N^ — 1) / (2 Nc), and n/ is the number of active flavors. 

Note that the basis in (9) is infinite, but the sums in (11) run, at a definite 
order in the loop expansion, over a limited set of operators. For instance at one 
loop the sub-matrix 7 of operators Oi/2,n results as follows, omitting a factor 
a,/(47r) 




(l,n) 
(2,n) 



(1, n-2) 

(6— n) (n— 5) 

(n-6) (n-5) 
2 

(2, n 
(n — 6) (n — 5) 

(4Cf-CA)(n-6)(n-5) 
2 



(1, n) 

4CF((l-n)(n-3)) 
3(-2+4?i-?i2) 

2) 



(1, n + 2) 

-((l+n)(2+n)) 

Ca 

(1+n) (2+n) 
2 



(2,n) 


(3Ca-4Cf) (2-4n+?i2A4c^ 



(2, n + 2) 
(1 + n) (2 + n) 

(4CF-CA)(l+n)(2+n) 
2 



(13) 



That means that even if the basis contains only operators 0(^x,i), 0(x,3) at the 
Mw scale, the evanescent operators 0(^^5) appear at one- loop order. 

The last step is to project the RG equations in 4 dimensions. This is most 
easily done by using the reduction formulas [24] to reexpress the insertion of 
evanescent operators in Green functions as a contribution to the coefficients of 
relevant operators 



Nm=J2r,,,N[R, 



(14) 



where the f matrix can be determined perturbatively. One is then able to decou- 
ple the evanescent operators at the level of the RG equation and this is equivalent. 



6 



as one can easily check, to define a non-minimal subtraction scheme which sets 
to zero the matrix elements of evanescent operators, as in [11]. The anomalous 
dimension matrix for relevant operators is modified as follows (schematically) 



7r 



7r,r' + 



/ , '~ir, e' e,r' 
e 



^J p. 



(15) 



At leading logarithmic order, only the reduction over magnetic momentum oper- 
ators is relevant, resulting in the following entries 



8 .. (3,5) 



X 



5 
1 



6 

r 



15 (4, 5)\Ca 



(16) 



the reduction over four fermion operators is of order a^ and it is relevant only for 
the NLO computation. 

We give the final results of the anomalous dimension matrix after the reduction 
to 4 dimensions in the basis (10), in a form suitable for comparison with Ciuchini 
et al. [11,12]: 



7 



^ I 7// 7/m 

47r V 7mm 



(17) 



where 7j/, ^mm are the scheme independent anomalous dimension matrices in 
the four-fermion and in the magnetic momentum sectors respectively. 



7// 



1 /- — 

' Ca 



V 



6 









Jn 



7 

7 (8Cf 

8 \8Cf 



2 
6 

6 
'Ca 








3 





3Ca 
22 



3Ca 
2(9CA-nf) 

3Ca 



2"/ 



3Ca 





16Cp-4C/ 



4 


2 
3 

22 
3 
2(9-nfCA] 
3Ca 



3 



5 




2 


3Ca 

4 


3Ca 

2nj 


3Ca 
6 

Ca 

2n, 



3Ca 



6 




2 
3 
4 
3 
2n 
~3 



/ 



-6 

2(n/-18CF) 



:ii 



while 7jm is the scheme dependent matrix connecting dimension 6 and dimension 
5 operators 



7/^ 



1 / 
2 

3 

4 
5 
6 



32 



3Ca 
+ 



32 
3Ca 
40 Ca Cf 





232 Cf 

9 
32 Cj^ I 88 Cp 

3 "•" 9 

32 Ca _ 40 Cf 

3 3 

4Cf 
9 



-8Ca _|_ 92 Cf 
3 "•" 9 



+ GUf 



32 I 32 Cj^ I 88 Cp 
■3Ca ~'~ 3 "•" 9 
-4 - 8|a ^^ + 92|V ^^ 

32 _ 40 Cf 
3 



(4n/ + 27%) 



3Ca 



9 

- Guf 

124 Cf^ 
9 ""/ 



(19) 



obtained from the computation of the graphs in Fig. (1) for different insertions 
of 4 fermion operators. 

Following the authors of [11] we define nf = u + d, nf = d — 2u and u, d are 
the number of active up and down flavors. 

As shown in [11], the matrix 7/^ is scheme dependent, and in order to give a 
sensible result we have to compute the matrix elements of operators Q^^4^^q- 




Qt 



sg 



sg 



(Ql\ 
[QlJ 


(Ql\ 




(20) 



We have checked that the physical result is indeed scheme independent. The 
easiest way to show this in our context is to perform a finite renormalization 



N' [Qt 



F] N 
'i-j 



Q- 



(21) 



which sets to zero the matrix elements in (20). Note that this finite subtraction 
is independent on the coupling as, a consequence of the fact that the scheme 
dependence stems from penguin diagrams at zeroth order in QCD [11]. 
It is well known that the ADM matrix is modified as follows 



7 



F^F 



-1 



(22) 



Using the results in (20) to define the matrix F, the scheme dependent sub-matrix 
becomes: 





D 





B 










K 




FIG. 1. Two loop graphs needed for 



Ifn 





FIG. 2. Example of graphs related by generalized Fierz identities 



imf 



1 / 

2 
3 

4 
5 
6 



7 


208 Cf 

9 
232 Cf 

-16 Cp 
^^Uf- 12 CpUf 



12 



\- 



8 
6 

11^_4C^ 

ACA-lQCp-Quf 



(23) 



%; 



which as expected coincides with the result in the HV scheme [11], where the 
matrix elements in (20) are zero. 

IV. CONCLUSIONS 

We have confirmed the results of Ciuchini et al. [11,12] by using a method 
which appears well suited for the NLO order computation. 

We stress that our considerations are of practical nature: in the evaluation 
of QCD corrections to electroweak processes this technique appears less involved 
than the t'Hooft-Veltman scheme, without introducing any potential ambiguity 
in the regularization. 

Moreover this method presents other advantages, which will be elucidated 
in an extended description of the computation [21]: we just mention here that, 
thanks to generalized Fierz identities [22], graphs like the ones in Fig. (2) are 
related and the knowledge of the "bare" graph with an open loop (on the left 
side of Fig. (2)) for arbitrary values of n allows to determine directly the value 
of the traced graph. 

This and similar considerations will be very useful in the NLO computation, 
certainly a formidable task, but in our opinion worthwhile in order to test as 
stringently as possible the Standard Model predictions. 

ACKNOWLEDGMENTS 

One of us (G. R) wants to thank Prof. H. Georgi for interesting discussions. 



10 



REFERENCES 



Electronic address 
Electronic address 
Electronic address 
Electronic address 



cella@sunlO.difi.unipi.it 
curci@nivxpil.difi.unipi.it 
ricciard@physics.harvard.edu 
vicere@sunlO.difi.unipi.it 
[1] Ammar et al. , Phys. Rev. Lett. 71 (1993) 674. 
[2] Peter Kim, talk given at McGill Heavy Flavour Conference 1993, 

CLNS/93/1254. 
[3] B. Grinstein, R. Springer and M. B. Wise, Phys. Lett. B 202 (1988) 138. 
[4] R. Grigjanis, P. J. O'Donnel and M. Sutherland, Phys. Lett. B 213 (1988) 

355. 
[5] R. Grigjanis, P. J. O'Donnel and M. Sutherland, Phys. Lett. B 224 (1989) 

209. 
[6] B. Grinstein, R. Springer and M. B. Wise, Nucl. Phys. B 339 (1990) 269. 
[7] P. Cho and B. Grinstein, Nucl. Phys. B 365 (1991) 279. 
[8] G. Cella. G. Curci, G. Ricciardi and A. Vicere, Phys. Lett. B 248 (1990) 181. 
[9] M. Misiak, Phys. Lett. B 269 (1991) 161. 
[10] M. Misiak, Nucl. Phys. B 393 (1993) 23. 
[11] M. Ciuchini, E. Franco, G. Martinelli, L. Reina and L. Silvestrini, Phys. Lett. 

B 316 (1993) 127. 
[12] M. Ciuchini, E. Franco, L. Reina and L. Silvestrini, preprint ROME 93/973 

(1993). 
[13] M. Misiak, preprint TUM-T31-46/93 (1993). 
[14] A. J. Buras, M. Misiak, M. Miinz and S. Pokorski, preprint TUM-T31-50/93 

(1993). 
[15] G.'t Hooft and M. Veltman, Nucl. Phys. B 44 (1972) 189. 
[16] W. Siegel, Phys. Lett. B 48 (1993) 193. 
[17] A. J. Buras, M. Jamin, M. E. Lauthenbacher and P. H. Weisz, Nucl. Phys B 

370 (1992) 69. addendum-ibid. B 375 (1992) 501. 
[18] M. Ciuchini, E. Franco, G. Martinelh and L. Reina, Phys. Lett. B 301 (1993) 

263. 
[19] L. V. Avdeev, G. A. Chochia and A. A. Vladimirov , Phys. Lett. B 105 

(1981) 272. 
[20] G. Curci and G. Ricciardi, Phys. Rev. D 47 (1993) 2991. 
[21] G. Cclla, G. Curci, G. Ricciardi and A. Vicere, in preparation. 
[22] L. V. Avdeev, Theoreticheskaya i Matematicheskaya Fizika 58 (1984) 308. 
[23] A. D. Kennedy, J. Math. Phys. 22 (1981) 7. 
[24] G. Bonneau, Nucl. Phys. B 167 (1980) 261, 
Nucl. Phys. B 171 (1980) 447. 



11 



